Analytical model of the response of a superconducting film to line currents 
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We theoretically investigate the response of a superconducting film to line currents flowing in 
linear wires placed above the film, and we present analytic expressions for the magnetic-field and 
current distributions based on the critical state model. The behavior of the superconducting film 
is characterized by the sheet-current density K^, whose magnitude cannot exceed the critical value 
jcd, where jc is the critical current density and d is the thickness of the film. When the transport 
current 7o fiowing in the wire is small enough, \Kz\ is smaller than jcd and the magnetic field is 
shielded below the film. When Jo exceeds a threshold value 7co oc j^d, on the other hand, \Kz\ 
reaches jcd and the magnetic field penetrates below the film. We also calculate the ac response of 
the film when an ac transport current flows in the linear wires. 

PACS numbers: 74.25.Sv, 74.25.Nf, 74.78.-w 



I. INTRODUCTION 

The response of superconducting films to homogeneous 
applied magnetic fields is well understood, and analytic 
expressions for the distributions of the magnetic field 
and current density have been derivediiSiSiii^ based 
on Bean's critical state model.^ When small current- 
carrying coils are placed near the surface to probe the 
local properties of superconducting films, the magnetic 
fields generated by the coils are inhomogeneous. An- 
alytic expressions describing the response of supercon- 
ducting films to small coils have been derived for the lin- 
ear response regime j^i^iifliiLiSii^ but in order to measure 
the local distribution of the critical current density jc in 
superconducting films, it is necessary to investigate the 
nonlinear response ji^ti^iiSiiLiS Numerical computations 
of the nonlinear response of superconducting films to the 
inhomogeneous magnetic fields arising from small coils 
have been carried out in Refs. .11,19 20, 21 22i. but here 
we present analytic results for the nonlinear response to 
line currents above superconducting films. 

The procedure proposed by Claassen et al}^ for in- 
ductive measurements of the local jc distribution in films 
of thickness much greater than the London penetration 
depth A can be described briefly as follows. A small 
coil carrying a sinusoidal drive current /q cos ut is placed 
just above a superconducting film, and the induced volt- 
age V{t) = Vn cos{nLot + -dn) in the coil is detected. 
The amplitude of the third-harmonic voltage V3 is mea- 
sured as a function of the drive current amplitude Iq, 
and the threshold current Icq is defined such that V3 = 
for < Iq < Ico and V3 > for Iq > IcO- Because 
IcO oc jcd, where d the film thickness, jc can be evaluated 
from Icom^ii^^ Since an electric- field criterion must be 
applied for a precise determination of jc, it is important 
to evaluate the electric field Ef generated in the super- 
conducting filmiS 

In the present paper we consider linear wires as simple 



models of coil wires, and we analytically investigate the 
response of a superconducting film to linear wires carry- 
ing transport currents. In Sec. [H] we investigate the dc 
(ac) response of a superconducting film to a linear wire 
carrying a dc (ac) transport current: we determine the 
threshold current Ico, and we present the voltage V{t) 
and the harmonic voltages induced in the linear wire, as 
well as the electric field Ef induced in the superconduct- 
ing film. In Sec. IIIII we consider a superconducting film 
and two linear wires carrying transport currents of op- 
posite directions. We briefly summarize our results in 
Sec. El 



II. A SINGLE LINEAR WIRE AND A 
SUPERCONDUCTING FILM 

In this section we consider a superconducting fllm and 
a linear wire carrying a transport current, as shown in 
Fig. An infinitely long wire, parallel to the z axis, is 
situated at {x,y) = (0,j/o) where j/o > 0. The radius of 
the wire is assumed to be much smaller than yo. A 
superconducting film, infinitely extended in the xz plane, 
is situated at ~d/2 < y < +d/2, where the film thick- 
ness d is much smaller than y^ but is larger than the 



FIG. 1: Configuration of a superconducting film at y = and 
a linear wire at {x, y) = (0, yo)- The film is infinitely extended 
in the xz plane, and the infinite wire is parallel to the z axis. 
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London penetration depth A. Flux pinning in the film is 
characterized by the critical current density jcj which is 
assumed to be constant (independent of magnetic field) 
as in Bean's critical state model^ and to be spatially ho- 
mogeneous in the film. We consider the limit d ^ 0, 
as this simplification allows us to obtain simple analytic 
expressions for the magnetic-field and current distribu- 
tions. In the thin-film limit of d — > 0, the sheet current 
Kz{x) — J^^l^ jz{x,y)dy plays crucial roles, and the up- 
per limit of \Kz\ is the critical sheet-current density jcd. 
The lower critical field Hci is assumed to be much smaller 
than jcd (i.e., Hd/ jcd 0), such that the flux penetra- 
tion into superconducting films is not affected by -ffci, 
but is determined by jcdP''^ 

We introduce the complex field 7i(C) = Hy{x,y) + 
iHx{x,y)r^'^^^ which is an analytic function of C = 
x + iy for y ^ and (x, y) ^ (0, yo). The Biot-Savart law 
for the complex field is given by 



1 /■+°° 
^o(c) + ^y du 



(1) 



where 7io(C) is the complex field arising from the line 
current alone. The 7io(C) is given by 



Ho(C) 



27r C - iyo ' 



(2) 



where Iz is the transport current flowing in the linear 
wire. At the upper ((^ = x + ie) and lower ~ x ~ ie) 
surfaces of the superconducting film, where e — d/2 is 
a positive infinitesimal, the perpendicular and parallel 
magnetic fields Hy{x,0) = ReH{x±ie) and Hrc{x,±€) — 
lmH{x ± ie) are obtained from Eq. as 



Hy{x,0) 



ReHoix) + — P 

ZTT 



H,{x,±e) = linno{x)TKz{x)/2, 



du^, (3) 

X — u 

(4) 



where P denotes the Cauchy principal value integral. The 
complex potential is defined by Q(() = J H{()d(, and 
the contour lines of the real part of G{C) correspond to 
magnetic-field lines. 

The magnetic flux per unit length around the linear 
wire is 



-Mo / dyH^{0,y) 

'ya+ro 



/ioRe 



e(«(yo + ro)) + lim g{iv) . (5) 



We have introduced a cutoff length tq ^ yo, where rg 
is of the order of the radius of the wire, to remove the 
logarithmic divergence as C ~* iyo- The magnetic flux 
per unit length <I>j(x') up through the film {y = 0) in the 
region x' < x < +oo is 

/"OC 

$/(a:') ^ Ho duHy{u,0) 



- g{x') + lim g{u) 



(6) 



A. DC response 

In this subsection we consider the magnetic-field dis- 
tribution when the linear wire carries a dc current Iz = 
/q > that is held constant after monotonically increas- 
ing from Iz — 0. 



1. Linear response for < /o < /cO 

For < /q < /cO, the magnetic field is completely 
shielded below the film, y = ImC < 0. The field distribu- 
tion can be obtained by the mirror-image technique, and 
the resulting complex field is 



.^^-^ forImC>0, 
for Im C < 0. 



(7) 



The complex potential Q{C) — J 7Y(C)<iC for ImC > is 
given by 



g{C) = — iarctan ( — 
71" V^o 



(8) 



The perpendicular magnetic field and sheet-current den- 
sity are thus given by Hy{x, 0) = and 



Kz{x) - - 



Iq yo 

9 1 9 ' 



(9) 



respectively. The net current induced in the supercon- 
ducting film is Kz{x)dx — — /q, as expected. Note 
that the sheet-current density cannot exceed jcd; that is, 
1-^2(2^)1 < jcd. Because the maximum of \Kz\ given by 
Eq. Q is /o/ttj/o, Eq- ® is valid for Iq < Ico, where the 
threshold current is given by 



IcO = TTjcdyo- 



(10) 



Figure 121 a) shows the magnetic- field lines [i.e., the con- 
tour lines of Re5(a; -I- iy)] calculated from Eq. (|SJ, and 
the dashed line in Fig.|2Id) shows Kz{x) given by Eq.©. 

The magnetic flux per unit length around the linear 
wire, calculated by substituting Eq. {Tjl into Eq. is 



dy ■ 



yo 



yo+ro 



y^ 



— 2 ^ ^ola, 
Vo 



(11) 



where the inductance per unit length Lq = 
(^o/27r) ln(27/o/''o) corresponds to the difference be- 
tween the self inductance of the linear wire and the 
mutual inductance of the wire and its image. Because 
the perpendicular magnetic field in the film is zero, the 
magnetic flux up through the film defined by Eq. JH)) is 
also zerotSS 



^f{x) = 0. 



(12) 



Equations IJI), ijSJl, ifTTIl . and ifT^ are vahd for 
< /o < ho- 
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FIG. 2: (color online) Magnetic-field and sheet-current distributions (a,d) for Iq/Ico = 0.99, (b,e) for Iq/Ico = 2 (a = 1.73), and 
(c,f) for Iq/IcQ — 3 {a — 2.83). Top figures (a,b,c) show the magnetic-field lines around a linear wire and a superconducting 
film. The cross symbols at (a;, y) = (0, j/o) ~ (0, 1) denote the position of the linear wire, the black horizontal line the shielded 
region > a of the superconducting film, and the gray horizontal lines the penetrated region \x\ < a of the film. The 
vertical dot-dashed lines indicate the positions of the flux front at x = ±a. Bottom figures (d,e,f) show the distributions of the 
perpendicular magnetic field Hy{x,0) divided by jcd (solid lines) and the sheet-current density Kz{x) divided by jcd (dashed 
fines) in the superconducting film. 



2. Nonlinear response for Iq > Ico 



IS 



For Iq > Jco, on the other hand, the maximum \Kz\ 
reaches jcd and the magnetic field penetrates below the 
superconducting film. The field distributions for Iq > 
Ico, therefore, must satisfy 

Kz{x) = -jcd and Hy{x,0) ^ for < a, (13) 
\Kz{x)\ < jcd and Hy{x,0) = for |a;| > a, (14) 

where the flux fronts [i.e., the boundaries between the 
region of Hy{x,0) ^ and that of Hy{x,0) = 0] are at 
X — ±a. The complex field and the parameter a are 
determined such that they are consistent with Eqs. H13() 
and H14|) . as derived in Appendix^ The parameter a is 
a function of /q, 



and the complex field is 



(15) 



where the upper sign holds for Im(^ > and the lower 
sign for Im^ < 0. The corresponding complex potential 



i arctan ( — I — arctanh 

yo, 



(17) 



Figures [2Ib) and (c) show the magnetic-field lines, i.e., 
the contour lines oi'ReQ{x + iy) calculated from Eq. I|17|l . 

The perpendicular magnetic field Hy{x, 0) = ReTi.{x± 
ie) in the superconducting film, obtained from Eq. H16|l . 
is 



' jcd xV a? - x^ 
i/,(x,0) = <; 2 x^^yl ^°'-|^l<«' (18) 
^0 for \x\ > a, 

and the sheet-current density Kz{x) = \m.\H.{x — ie) — 
H{x + ie)] is 



-jcd 



for \x\ < a, 



K^ix) 



-Jcd 1 - ^ — 1 for \x\ > a. 



(19) 



The net current induced in the superconducting film is 
again Kz{x)dx = —Iq, as expected. Figures EJe) 
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and (f) show Hy{x, 0) and Kz{x) given by Eqs. (|18() and 



B. AC response 

In this subsection we consider the time-dependent field 
distributions when the hncar wire carries a sinusoidal ac 
drive current Iz{t) = Iq cos cut. In inductive measure- 
ments of the local jc in superconducting films, harmonic 
voltages induced in coils are detected li^^iSii^ For precise 
jc measurements, it also is important to determine the 
electric field induced in the filmi^ 

We wish to calculate the magnetic field around the 
linear wire defined by Eq. and the voltage 

per unit length induced in the wire, 

V{t) = R^Iq cosujt - d^.^{t)/dt, (24) 

where Rw is the resistance per unit length of the wire. 
We also wish to determine the magnetic flux per unit 
length up through the film defined by Eq. lO, 

and the electric field in the filmjSS 

Ef{x, t) = -d^f{x, t)/dt. (25) 

1. Linear response for < /o < /cO 

For < /o < IcQ the magnetic field is completely 
shielded below the superconducting film, y = Im^ < 0, 
as in Sec. Ill A II The complex field, the sheet-current 



(fT^ . respectively. 



(20) 



(21) 



(22) 



(23) 



I 

density in the film, and the magnetic flux around the 
linear wire are given by Eqs. {Tj), and (|ll|l . respec- 
tively, except that now /q in those equations is replaced 
by /q cos ujt. 

The magnetic flux per unit length around the linear 
wire is given by $u,(t) = LqIq cos ujt. The voltage in- 
duced in the wire defined by Eq. (|24|l is thus given by 
V{t) = RwIoCOSLut + ujLoIosinujt. For < /q < IcO, 
the harmonic voltages, the magnetic fiux per unit length 
penetrating the film $ / , and the electric field in the film 
Ef are all zerom^ 



2. Nonlinear response for Iq > IcO 

For Iq > /co, on the other hand, the magnetic field pen- 
etrates through the superconducting film, as discussed in 
Sec. lirX2l 

For ac drive current Iq cos tot, the magnetic flux per 
unit length around the linear wire <I>^ (t) is'* 

$»(0 = $»(/o) - 2$„(/o(l - cosLut)/2^ (26) 

for <ujt <TT, and ^w{t) = -^w{t- tt/uj) for n < ut < 
2tt, where $u,(/o) is defined by Eq. (|21|l . The voltage per 
unit length of the wire, calculated from Eqs. H24|l and 



The magnetic flux per unit length around the linear wire, calculated by substituting Eq. (|17|l into Eq. jSJ, is 



^1, 



a^ + yo + yo? + Vo In 



2/6 



Combining Eqs. (dJ, and yields $^ = *to(-?^o)7 where 



r r 



IcO - /o + -^0 In ( 4^ 



'cO 



for 0< Iq< IcO, 
for Iq > IcO- 



The magnetic flux up through the fllm, calculated by substituting Eq. (|17|1 into Eq. @, is 



<^f{x) 



■ \/ o? — x^ + ^jo? arctanh | W ^ 



for |a;| < a, and $/(a;) = for |a;| > a. Combining Eqs. (^^J, CSll, and yields ^f(x) = $/(a;,/o), where 
for < /q < /cO or |a;| > a. 



$/(x,/o) = 



Mo 
27r 



/ /q - ( 1 + ) + /oarccosh 



a;2 + J/2 



for Iq > /co and \x\ < a. 
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(ESli, 



Vit) 



RwIq cos cut 



(/o(l -costji)/2), (27) 



where L^{Iq) = d$u,(/o)/d/o is the differential induc- 
tance given by 




for < /o < /co, 
for lo > I co- 



rn 



In response to the ac drive current, the magnetic flux 
per unit length up through the film $y(a;,i) ia^ 

= $/(x,/o) - 2^f(x,lQ{l - cosujt)/2^ (29) 

for < ujt < TT, and $y(a;,t) = ~^f{x,t - tt/lo) for 

IT < ujt < 2tt, where $/(x,/o) is defined by Eq. (j^ . 
The electric field induced in the film, calculated from 
Eqs. (Eni and is 

Ef{x,t) = w/o sino^iL/^a;, /o(l — cosa;i)/2^, (30) 
where the function Lf{x,Io) = f{x, Io)/dIo is given 

by 



for < /o < IcO or > a, 

yo \ (31) 



i/(^,/o)= ^^arccosh A ^ 

for /q > IcQ and |x| < a. 

In order to measure jc in superconducting films by the 
inductive method detecting the harmonic voltageS(i^ it 
is important to estimate the induced electric field in su- 
perconducting filmSfS^ The maximum electric field is in- 
duced just below the linear wire at a; = 0, and is given 
approximately by 

for < /o - /cO < IcO- 

3. Harmonic voltages 

The voltage per unit length induced in the linear wire 
can be expressed as the Fourier series 



V{t) =J2VnCOs{nUjt + dn), 



(33) 



n=l 



where the amplitude Vn and phase difference of the 
nth harmonics are calculated by 

K exp(-ii?„) = - / d{u;t)V{t)exp{inujt). (34) 



0.3 



g 0.2 
_^ 

> 0.1 



- 


1 1 1 1 


/ / 
* / 

j/^ 









0" 



-90° 
-180" 



1 2 3 4 5 

FIG. 3: (color online) The amplitude V3 and the phase -dz of 
the third-harmonic voltage as functions of the drive current 
amplitude Jq. 



Because of the periodicity of V(t -\- tt/lu) = —V{t), the 
even harmonics are zero; i.e., Vn — Q for n = 2, 4, 6, . . .. 
The fundamental voltage 14 exp(— i'i?i) is simply deter- 
mined by the resistance per unit length of the linear wire, 
the self-inductance per unit length of the wire, and the 
mutual inductance per unit length between the wire and 
its image. We are interested in the odd harmonics of 
Vn exp(— z-iJn) for n = 3, 5, 7, . . .; they are obtained by 
substituting Eqs. (|2Z|) and (|2Hll into Eq. which yields 

I4exp(-ii3„) 



— cj/q / (i6'exp(m6') sin^L^.I /o(l — cos6')/2 



^w/o / exp (in6') sin6' In [—^■ 



/o 1 — cos t 



,(35) 



where 9c = arccos(l — 2Jco//o)- ForO < /q/Ico — 1 ^ 1 
we have 



Vn 



fJ-O 



bjl. 



IcO 

16n / /o 



1 



1/2 



(36) 
(37) 



Figure O shows the /q dependence of the third- 
harmonic voltage V3 exp(— 11^3) calculated from Eq. H35|l 
with n = 3. Although the present model of a linear wire 
is oversimplified, the behavior of the third harmonic volt- 
age shown in this figure qualitatively agrees with the ex- 
perimental data measured by a coil with a YBa2Cu307_y 
film,i^*iLia 



III. TWO LINEAR WIRES AND A 
SUPERCONDUCTING FILM 

Two linear wires are better than a single wire to model 
current-carrying coils. In this section we consider a su- 
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perconducting film and two linear wires carrying trans- 
port currents of opposite directions, as shown in Fig. 0] 
Infinitely long wires, parallel to the z axis, are situated 
at (x^y) — {±xo,yo) where xo > and yo > 0. A su- 
perconducting film, infinitely extended in the xz plane, 
is situated at ~d/2 < y < +d/2, where A < d <C yo- 

One wire carrying a dc current +1^ is situated at 
{x,y) = {+xo,yo) a-nd another wire carrying a current 
of opposite direction —1^ is at {x,y) — (— a;o,?;o), where 
xq > and j/o > 0. 

The complex field due to the two wires is given by 



Wo(C) 



2tt 



1 



c- 



where Ci = + iyo and C2 



Ci C - C2 
-xo + iyo = 



(38) 



A. DC response 



In this subsection we consider the magnetic-field dis- 
tribution when the linear wire carries a dc current 1^ — 
Iq > that is held constant after monotonically increas- 
ing from Iz — 0. 



1. Linear response for < /o < /cO 

For < /o < /co, the magnetic field is completely 
shielded below the film, y = Im^ < 0. The field distribu- 
tion can be obtained by the mirror-image technique, and 
the resulting complex field is 



27^11 (C) forImC>0, 







where 



W||(C) 



c 



for ImC < 0, 



c 



The complex potential for Im {Q > is given by 



(39) 



(40) 



(41) 









t) 


--Vo 


+-Vo 1 



FIG. 4; Configuration of a superconducting film at y = 
and two linear wires at {x,y) = {±xo,yo) carrying transport 
currents ±7^. 



The perpendicular magnetic field is Hy{x,Q) = 0, and 
the sheet-current density is thus given by 



Kz{x) = -loF^ix) 



(42) 



where the function Fw{x) is determined by the configu- 
ration of the wires, 



F^{x) 



yo 



yo 



{x~xo)'^+yo {x + xoy+yo 



(43) 



The sheet-current density [^'^(x)! is maximum at a; = Oc, 
where 

The sheet-current density cannot exceed jcd; that is, 
|-f'^z(2^)| < IoFw{ac) < jcd, and Eq. is valid for 

-^0 5: IcO, where the threshold current is given by 



IcO = jcd/F^{ac). 



(45) 



Figure Ela) shows the magnetic-field lines [i.e., the 
contour lines of KeG{x + iy)] calculated from Eq. H41|l . 
and the dashed line in Fig. |3[d) shows Kz{x) given by 
Eq. 

2. Nonlinear response for Iq > Ico 

For Iq > Ico, on the other hand, the maximum \Kz\ 
reaches jcd and the magnetic field penetrates below the 
superconducting film. The field distributions for Iq > 
Ico, therefore, must satisfy 

Kz{x) — —jcdsgn{x) for b < \x\ < a, (46) 
Hy{x,0) = for l^l < 6 or > a, (47) 

where the fiux fronts are at a; = ±a and x — ±b. The 
complex field and the parameters a and b (where a > 
6 > 0) are determined such that they are consistent with 
Eqs. (gnil and ^7^, as derived in Appendix EJ 

The parameters for flux fronts, a and 6, are determined 
as functions of Iq for Iq > /co, by solving the following 
two equations. 



aE(k) K(k), (48) 

a 




a[K{k) - E{k)], 



(49) 



where K{k) and E{k) are the complete elliptic inte- 
grals of the first and second kinds, and k = ^1 — b"^ /a?. 
These parameters obey a — b — ac when Iq = Icq, and 
a > ttc > b > Q when Iq > Icq, where ac is given by 
Eq. (|44|) . As shown in Fig.|Sl a increases and b decreases 
with increasing Iq. 
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(a) 








_/o//co=0.99 


(d) 


\ 

\ 

\ 


//v(A-,0) 




-/o//,„=0.99 





-2 



+2 +4 -4 




FIG. 5: (color online) Magnetic-field and sheet-current distributions (a,d) for lo/ico = 0.99 (ttc = 1.07), (b,e) for lo/ico ~ 1.5 
(a = 2.06 and b — 0.32), and (c,f) for lo/lao = 3 {a — 3.23 and b = 0.02). Top figures (a,b,c) show the magnetic-field lines 
around linear wires and a superconducting film. The cross symbols at {x,y) = (±a;o,j/o) = (±1,1) denote the positions of 
the linear wires, the black horizontal lines the shielded regions (|a;| < fo or \x\ > a) of the superconducting film, and the gray 
horizontal lines the penetrated regions & < < a of the film. The vertical dot-dashed lines indicate the positions of the flux 
fronts at X = ±a and ±6. Bottom figures (d,e,f) show the distributions of the perpendicular magnetic field Hy{x,0) divided by 
jcd (solid lines) and the sheet-current density Kz{x) divided by jcd (dashed lines) in the superconducting film. 



The complex field is derived in Appendix [HI and is 
jiven by 



(50) 




(A-o,yo)=( 1,0.3) 



1 1.2 1.4 1.6 1.8 2 

/0//c0 

FIG. 6: (color online) The parameters for flux fronts, a (upper 
lines) and b (lower lines), as functions of Jo, determined by 
Eqs. 14811 and 1491 . The solid lines show a and b for {xo,yo) ~ 
(1,3) (oc = 1.92), the dashed lines for {xo,yo) = (1, 1) (nc = 
1.07), and the chained lines for (a;o, yo) = (1, 0.3) (bc = 1.00). 



where 
W±(C) 



1 



1 



du 



(51) 
(52) 

(53) 



Integrating Eq. H5U|) [or ljC9p ]. we obtain the complex 
potential CJ(C) — jH{C)d(^, as 



An 



e-c! 



+ ^[-Go(C,Ci) + Go(C,C2)] 



— / duGoiCu) 

TT Jb 



where 



Go{C,s) = arctanhi 



q2 „ (^2 g2 _ ^2 



(54) 



(55) 



Figure|S{b) and (c) show the magnetic-field lines, i.e., 
the contour lines oiReG{x + iy) calculated from Eq. ltM|l . 
Figure |S{e) and (f) show Hy{x,0) = ReH{x ± ie) and 
Kz{x) = Im['H{x — ie) — T-l{x + ie)] calculated from 
Eq. ^ [or from Eqs. g^), E3, Ip^ . and fHTO^ ]. 
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B. AC response 

The response of a superconducting film to ac drive cur- 
rents zklz (t) — ±Iq cos U!t flowing in the two linear wires 
at {x,y) = (±xo,yo) is calculated in a manner similar to 
that described in Sec. IIIBI The magnetic flux per unit 
length <i>u, linked in the two wires is calculated from 

/+xo—ro 
dxHy{x,yo) 



= /ioRe GiC2+ra) -GiCi-ro) 



(56) 



and the magnetic flux per unit length ^f{x') up through 
the film (y = 0) in the region x' < x < +oo is calculated 
from Eq. ©, where Q{() is given by Eqs. ||1T)) and 
The time-dependent magnetic flux is given by Eqs. H2ti|) 
and The voltage V{t) induced in the wires per unit 
length and the electric field induced in the film are given 
by Eqs. and (|2SJ), respectively. 



Eqs. 0, dl, (EH), and |(in|)- The behavior of the film 
can be summarized as follows. 

(i) Response to a dc current, Iz — la'- For a small dc 
current in the range < Iq < Ico oc jcd, the sheet-current 
density in the film is less than the critical value jcd, 
and no magnetic field penetrates into the region y < 
below the film. However, for a large dc current Iq > 
/co, the sheet-current density reaches \Kz\ = jcd and the 
magnetic field penetrates into the region below the film. 

(u) Response to an ac current, 1^ = locosojt: For 
< Iq < Ico the voltage V{t) induced in the linear wire 
has only the first harmonic at the fundamental frequency 
LUi — Lo. However, for /q > /co, odd-harmonic voltages 
with frequency w„ — nuj, n = 3, 5, 7,..., are also induced 
in the linear wire. 

The complex field for the case of a single current- 
carrying linear wire above the film is derived in Appendix 
A. The complex field for the case of a pair of current- 
carrying wires in a plane perpendicular (parallel) to the 
film is given in Appendix B (Appendix C). 



IV. DISCUSSION 

Here we discuss the more realistic situation of exper- 
iments detecting the response of a superconducting film 
to a circular current-carrying coilJ^ We consider a single- 
turn coil parallel to a superconducting film. When the 
radius of the coil is xq and the distance between the coil 
and the film is yo, the configuration of the coil and the 
film in the xy plane is similar to that in Fig. ^ We 
expect that the magnetic-field lines, the magnetic-field 
component perpendicular to the film, and the circular 
sheet-current distribution in the film then will be similar 
to the corresponding quantities shown in Fig. O 

In actual superconducting films the critical current 
density jc and the film thickness d can be inhomoge- 
neous, although in the present paper we assumed that 
jcd is spatially homogeneous. The magnetic field pro- 
duced by the coil of radius a:o is largest directly below 
the coil (i.e., the circular region of radius ~ xq in the 
film) when xq yo. If jcd in the circular region below 
the coil is nearly homogeneous, the response of the film 
will be similar to that presented in Sec. lIIII and the signal 
in the coil will yield the local value of jcd. In this sense 
the resolution of the measurements using the coil is on 
the order of xq. However, if jcd is very inhomogeneous 
within the circular region below the coil, the response of 
the film will be quite different from that given in Sec. lIIII 



V. CONCLUSION 

We investigated the response of an infinite supercon- 
ducting film at y = to linear wires above the film 
carrying transport currents 1^. We derived analytic ex- 
pressions for the complex field, and these are given in 
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APPENDIX A: COMPLEX FIELD FOR A 
SUPERCONDUCTING FILM WITH A LINEAR 
CURRENT-CARRYING WIRE ABOVE IT 

In this appendix we derive Eqs. (|15|l and (|16|l . Con- 
sider the function f{() defined by 



/(C)= [H(C)-r7||(C)] 
where 7]|| (C) is 

-^0 iyo 



vwiO 



^Tre+yl 



(Al) 



(A2) 



We calculate the integral /pdC'/(C')/(C' — C) along the 
closed contour C shown in Fig. |3 as 



C'-C 2ni J_ 



u-C 



(A3) 

where the contour integral along the infinite circle van- 
ishes because |/(C)| |W(C)CI ^ for |C| ^ oo. Be- 
cause the integrand in the left-hand side of Eq. ljA3|l has 
poles at C' = C and (' = iiyo, the contour integral of 
Eq. (|A3|) is also calculated by using the residue theorem, 
as /(C) - /o(C), and thus yielding 

fiO-foiO ^ TT- / du — , (A4) 



2711 



u-C 



9 




FIG. 7: (color online) The contour C in the i^' plane for the 
integrals in Eqs. 1A3II and IIU2t consists of a line just above the 
real axis, C,' = it + ie from u = —R' to u = +R where R' < R, 
an infinite circle, = 7?e'^ from 8 = e/R to 6* = 27r — e/R, 
a line just below the real axis, C^' = u — ie from u = +R to 
u = —R', and an infinitesimal line parallel to the imaginary 
axis, = —R' + iv from v — —e to v = +e, taking the limit 
of i? ^ <x, R' (X, and e — > +0. 



where /o(C) is given by 



/o(C) 



2n e + vl ■ 
For — X zt ie, the 7i(C) and ryu {() are given by 



(A5) 



H{x±ie) = Hy{x,0)+i[lmHo{x)TK,{x)/2], (A6) 
77||(a;) = ilmHo{x), (A7) 

respectively, and thus yielding 

n{x±ie) -r]\\{x) = Hy{x,0) T iK^{x)/2. (A8) 

By using Eqs. JJU, (|X8|) . and 



sgn(a;) vx^^^(? for |a;| > a. 



— x^ 



we have 

f{x + ie) - /(x - ie) 



for |a;| < a, 
(A9) 



— 2isgn(a;)-\/2;^ — a?Hy{x,Q) for jxl > a, 
a?- — a;^-ft'^(a;) for |a;| < a, 

'O for \x\ > a, 



a? — x^jcd for |x| < a, 



where we used Eqs. (|13|l and (|14|l . Substitution of 
Eq. iprUI) into Eq. (|X3I) yields 



/(C) -/o(C) 



jcd /■+° , Va^ - 
1^ — / du — 



(All) 



where the upper (lower) signs hold for Im C > (Im < 
0). 

From Eqs. IjAip . (|X5|) . and (|A11|) we obtain 



no 




(A12) 



Note that the first term of the right-hand side of 
Eq. (jA12p . which is proportional to Iq, corresponds to 
the complex field in the ideal Meissner state for two 
semi- infinite strips situated at |Re(^| > a with a line cur- 
rent at C = ii/Q. By using = T*C: we confirm 
that Eq. (|A12|I for a ^ is identical to Eq. Equa- 
tion (|A12|) can be rewritten as 



no 




(A13) 



(AlO) 



However, the last term of the right-hand side of Eq. (|A13|) 
must vanish, because 7i(C) is finite at C = ±a. We thus 
obtain Iq — njcdy^ + Uq, such that the parameter a de- 
scribing the position of the flux front is given by Eq. I|15|) . 
The resulting expression for the complex field is given by 
Eq. |TB1). 

See also Appendix IbI for the complex field for a super- 
conducting film with two wires in the yz plane. 



APPENDIX B: COMPLEX FIELD FOR TWO 
PARALLEL WIRES IN THE yz PLANE 

We present here the complex field for a superconduct- 
ing film with two parallel linear wires above it. An infi- 
nite superconducting film is in the plane y = 0, as shown 
in Fig.^ One wire carrying a dc current +Iq is situated 
at {x,y) = {0,yi) and another wire carrying current —Iq 
is at (x,y) = (0,2/2), where < yi < y2- The derivation 
of the complex field for the two wires is similar to that 
given in Appendix ^ for a single wire, and we exhibit 
here only the resulting expressions. 

For < /q < /co, the sheet-current density and the 
perpendicular magnetic field satisfy [-/^^(x)! < j^d and 
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Hy{x, 0) = 0, respectively. The threshold current is given 
by 



2/12/2 



2/2 - 2/1 

and the complex field for < /q < Ico is 



(Bl) 



for Im C > 0, 



for Im C < 0. 

(B2) 

For Jo > IcO, K^x) and Hy{x,0) fulfill Eqs. (O and 
(|14|l . and the relationship between a and /q is given by 



Jo = njcd 



(B3) 



The complex field and complex potential for /q > /co are 
respectively given by 



I zvr 



1 



1 



>2 , 2 I *2/i H , „ „ 

C'+2/?V + 2/? 



C2+2/I 



«2/2 



+ 2/2 



(B4) 



2 



2tt 



I arctan | — — arctanh i / ^ 

2/1 / \ V « + 2/1 



i arctan 



2/2 



arctanh 



/a2_C2 
+ 2/1 



(B5) 



APPENDIX C: COMPLEX FIELD FOR TWO 
PARALLEL WIRES IN A PLANE PARALLEL TO 
THE xz PLANE 

In this appendix we derive Eqs. (|48|l " (|52|l . Consider 
the function F{Q defined by 



^^(0= [H(C)-H||(C)] 0(C), 



(CI) 



where H||(C) and 0(C) are defined by Eqs. (|^ and 
respectively. We calculate the integral d('F{C')/{C' ~ 
C) along the closed contour C shown in Fig. [7| as in 
Eq. |J3^ . 



1 r,^,FiC) 



2m 



1 /■+°° , F{u + ie)- F(u-ie) 
du 



(C2) 

where the contour integral along the infinite circle van- 
ishes because |F(C)| - |7^(C)C^I for |C| -> oo. Be- 
cause the integrand on the left-hand side of Eq. ljC2p has 
poles at C' = C: iCi, ^iid ±C2, the contour integral of 



Eq. ljC2p can be calculated by using the residue theorem, 
as F(C) — -Po(C), ^iid thus yielding 



^^(C)-i^o(C) = ;^ 



The Fo(C) is defined by 



+ CXD 



du 



F{u + ie) - F{u - ie) 



(C3) 



MO - 



An 
27r 



0(Ci) 0(C2) 



-Ci) , 0(-C2) 



C-Ci C-C2 C + Ci C + C2 
C0(Ci) C0(C2 



(C4) 



where we used <?!'(— Ci) — — 0(Ci) and 0(— C2) — — 0(C2)- 
For C = a; ± «e, Eq. is reduced to 



ix ± ie) 



±z0(a;) 



for |a;| < 6 or LtI > a. 



[sgn(a:)0(a;) for h <\x\ < 



where 



0(x) = sgn(a — |x|)-\/|(a2 — x'^){x'^ — IP')] 

' ^J{a? ~ x'^)(lP - x"^) for < 6, 



(C5) 



(C6) 



■\/ [a? — x'^){x'^ — 6^) for 6 < < a, 
- \J (x^ — d?- ) (xp- — \P- ) for > a. 

Substituting Eqs. l|X8)l and (|C5|l into Eq. l(CT|l with C = 
X ± ze, we have 

F{x + ie) - Fix - ie) 

(2i(t>{x)Hy{x,0) for < b or |a;| > a, 

{—isgn{x)(l){x)Kz{x) for < |x| < a, 

{0 for |a;| < 6 or > a, 

ijcd(j){x) for & < |a;| < a, 

(C7) 

where we used Eqs. H46|l and (|47|l . Substitution of 
Eq. lEZj) into Eq. ^ yields 



nC)-^o(C) = ^/ 



0(u) 



j.d r,^.c0(.) 



such that 

W(C)=H||(C) 



^^o(c) , Jed c r^^. 



0(C) 0(C) A - C 



By using 

0(.) _0(C) , (c^-.^)(c^ + 



0(C) 0(s) 



0(C)0(S) 



(C8) 



(C9) 



(CIO) 
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we rewrite Eq. ljC9p as 

^(C)-H| 
= C<^(C) 



W(C)-H||(C) 

jcd 



1 



du 



c 



where 

no 



2ti 



no, 



(Cll) 



'^(Ci) 



jcd r , C' + - a' - 
/ du — 

Jb 4){u) 

- b^)n{±a) + - e)ni±b) 



62 



(C12) 



For C = ±a and ±6, we have 
7^(±a) 



27r 



2lT 





-62 












du\j 






/a2 





Jcd r 



du 







-62 






-CI 


u2 


-62 




a2 


— 






/a2 


-CI 






-62 


a2 


— u2 




u2 


-62 





(C13) 



(C14) 



In order to remove the divergences in Eq. (jClip at ^ = ±a 
and ±6, we require TZ{±a) — TZ{±b) — 0, thus yielding 



jcd 
jcd 



-Re^ 



Re^ 





-62 


a2 


-a 


'a2 


-a 



C2-62 



du 



du 



u^ 


-62 


a? 


9 

— 


a2 


— u2 


u2 


-62 



(C15) 
(C16) 



Equations HC15|) and (jC16p reduce to Eqs. ijlHl) and 
respectively. We thus obtain Tl{C,) = for any C, and 
Eq. l(nTT|) is reduced to Eqs. (EOJl, and 

The sheet current Kz{x) = lm[H{x — ie) — H{x + 
ie)] and the perpendicular magnetic field Hy{x,0) = 
Re7i(a; ± ie) are obtained by substituting C — x ± ie 
in Eq. ^ with Eq. l|J8)l . 



i7y (x, 0) T « i^2(a:)/2 = ± ie) - T^,, (x) 

X X 



ix±^e)^^^ 



jcd 



(x^^cmci) {x^-Q^)m) 

du 



TT Jfj [{x ± ie)"^ — u^]<f){u) 
x 



{x ± ie) I —Re 
jcd du 



(:^^-C?)'^(Ci] 

X 



- u) + S{x + u)^ 



(C17) 



Substitution of Eq. HC5|I into Eq. (|C17p yields the sheet 
current Kz{x) and the perpendicular magnetic field 
Hy{x,0). The Kz{x) for 6 < |x| < a is given by Eq. 
and that for < 6 or |x| > a is given by 



Kz{x) ~ —x4){x) I /qRc 

TT 



■jcd 



1 



du 



^ (a;2 — u'^)(/)[u) 



(C18) 



The Hy{x,0) for |a;| < 6 or |a;| > a is given by Eq. lH7|l . 
and that for 6 < |a;| < a is given by 



Hy{x,0) = -\x\(l){x) { loRe 

TT 



jcdP 



1 



{x^~eom) 

du 



/f, (x2 — U^)(j)(u) 

where P denotes the Cauchy principal value integral. 



(C19) 
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